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We study the normal forms for incompressible flows and maps in the neighborhood of an 
equilibrium or fixed point with a triple eigenvalue. We prove that when a divergence free 
vector field in R 3 has nilpotent linearization with maximal Jordan block then, to arbitrary 
degree, coordinates can be chosen so that the nonlinear terms occur as a single function of 
two variables in the third component. The analogue for volume-preserving diffeomorphisms 
gives an optimal normal form in which the truncation of the normal form at any degree gives 
an exactly volume-preserving map whose inverse is also polynomial inverse with the same 
degree. 

in 

1 Introduction 

j — . ' A system of ODEs £ = t>(£, t) gives rise to a volume-preserving flow when the vector field v has 

zero divergence, 

V-v = 0.L (f) 

• , 

Similarly, a diffeomorphism / : W l — > W l is volume preserving when its Jacobian has unit 
determinant, det Df(£) = 1. Volume-preserving dynamics arises, for example, in the flow of 



a Lagrangian element in an incompressible fluid [HolMl lDFG+861 IFKP881 ICFP961 ISAZMQ1L 



ISCVH041 IMJM051 IMJM06j . or the tracing of lines in a magnetic field, [TH851 ILF921 IGre93j . 
Thus such flows are of interest in the study of the motion of passive tracers in time-dependent, 
incompressible fluids and for the dynamics of charged particles in strong magnetic fields. Volume- 
preserving maps are a natural generalization of area-preserving maps to higher dimensions. They 
also arise as the normal form for certain homoclinic bifurcations for three-dimensional systems 
[(X)ST05llGMO06| . and as integrators for incompressible flows [QZ931 [Qm95l lKS95l ISnTM] and 
thus have intrinsic mathematical interest |CS90al IBKR901 IXia921 IRKKCA931 ILM981 IGM021 
ILM001ILM03] . 
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An equilibrium of a volume-preserving flow can undergo a number of codimension-one bifur- 
cations. For example when one of the eigenvalues vanishes and the remaining eigenvalues do not 
lie on the imaginary axis, then the resulting bifurcation is generically of standard saddle-node 
type. Similarly, when a single pair of eigenvalues lie on the imaginary axis, the bifurcation 
is of Hopf type. Since the eigenvalues of the equilibrium satisfy Y27=i ^ = there are two 
more exotic cases that are codimension-one only for low dimensions. For three dimensions, the 
codimension-one configuration {0, iu, —iuj} can give rise, through a Hopf-saddle node bifurcation 
to the creation of a periodic orbit surrounded by a family of invariant tori [Bro81]. The final 
codimension-one case occurs in four dimensions when the eigenvalues are — iu)2, w^2> —^2}', 
this bifurcation is analogous to the Hamiltonian-Hopf bifurcation [Bro81J. 

The linearization of a volume-preserving map at a fixed-point x* = f{x*) gives a matrix 
M = Df{x*) with determinant one, so that the multipliers satisfy Yii=i = For the three- 
dimensional case, the characteristic polynomial has the form 

p(A) = det(A/ - M) = A 3 - rA 2 + aX - 1 . (2) 

where r denotes the trace, r = trM, and a the second-trace, a = \ ((trM) 2 — trM 2 ). Con- 
sequently, (r, o") parametrize the space of volume-preserving matrices that are linearly conju- 
gate to M, at least when the eigenvalues are distinct. The eigenvalue configurations in this 
space, see Fig. Q] [LM98, ILMOOj . show that there are two codimension-two points. The first, 
at (t, a) = (—1, —1) where the eigenvalues are (1, —1, —1) corresponds to simultaneous period- 
doubling and saddle-node bifurcations, and the second, at (r, a) = (3, 3), to the triple eigenvalue 
(1, 1, 1). These are connected by a line segment {r = a : —1 < r < 3} on which two eigenvalues 
are on the unit circle, and the third is equal to one, corresponding to simultaneous saddle-node 
and Neimark-S acker bifurcations. 




Figure 1: Classification of the eigenvalues for a three-dimensional, volume-preserving map as a function of the 
trace r and second trace a. 

In this paper we study the normal forms for the codimension-two cases of a three-dimensional 
incompressble flow with eigenvalues {0, 0, 0} and of a volume-preserving map with multipliers 
{1, 1, 1}. We first recall some of the standard results on normal forms. 

Near an equilibrium, a set of ODEs with a smooth vector field takes the form 

£ = v(0 = jz + ko , (3) 
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for £ G M n , with a vector field u whose linear part is J£ and nonlinear part is &(£) = 0(2). A 
local normal form is a conjugate system 

?? = w(rj) = Jij + c(rj) , (4) 

that is "simpler" in some sense. For example, it is usually desirable to eliminate as many of the 
nonlinear terms in b as possible since then the vector field w will have fewer parameters. A first 
step in this process is to choose coordinates in which J itself is simple, and it is typical to begin 
by normalizing the linear part, that is to choose J to be in Jordan normal form. 

Since the linear parts of ([3]) and (j3|) are the same, the transformation can be assumed — to 
lowest order — to be the identity, 

v = i>(0 = t + HO . (5) 

Though the "simplest" form of a flow may only be topologically conjugate to the original system 
(as in the Hartman-Grobman theorem when J is hyperbolic), it is not possible to explicitly 
construct w unless we assume that both the vector fields and the transformation are smooth. 
When if) is a diffeomorphism then the vector fields of ([3]) and (JH) are related by 

Lhl>{Ov(£) = => 
b(C) + Dh(OJt + Dh(OKO) = JHO + c( JH + h(0) 

Under the assumption that v,w,h S C°°, the normal form can be computed by power series 
expansion. One way to accomplish this is to transform the terms of each degree one at a 
time. When this normalization has been carried out for all terms through degree d — 1, i.e., 
v = w + 0(d), we let h(£) = 0(d) in ©; these degree d terms must then solve the "homological 
equation" 

Cj(h)(0=c(0-K0 + O(d+1). (6) 
Here Cj is the homological operator defined by 

Cj = adj = [Ji, ■} , (7) 

where [, ] the Lie bracket of vector fields: 

adj h(0 = (J£ • V)h(0 - (h(0 ■ V) J£ = Dh{$JZ - Jh(0 . (8) 

The problem of constructing h at degree d then reduces to linear algebra on the finite dimensional 
space, F£, of n-dimensional vectors of polynomials in n- variables of homogeneous degree-d In 
this case, a solution h € to (|7j) exists providing c — b € rng£j n F2- If the operator C j 
were surjective on F£, then we could set c = to eliminate all nonlinear terms. However, the 
homological operator invariably has a nontrivial kernel, and c must be chosen to eliminate any 
terms in b that are not in rng£j. 

Recall that any matrix J = S + N has a unique splitting into commuting matrices such 
that S is semisimple and N is nilpotent. Moreover, ker J = kerS" n ker N. It can be similarly 
shown that Cj = Cs + Cn is a semisimple- nilpotent splitting of the homological operator on F^. 
The construction of a normal form requires the selection of a complement to rng£j; this can 
be chosen to be the intersection of complements to the ranges of the semisimple and nilpotent 
parts. 

For the semisimple case, it is easy to find a basis for F£ in which £$ is also diagonal (the 
vector monomials, see ^2]). In this case ker £5 is a complement to rngjCg, so that c can be 
chosen to be the terms in b that correspond to eigenfunctions of the homological operator with 
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zero eigenvalue. In this way the construction for the semisimple case reduces to characterizing 
ker £5. 

The construction of a complement to the range of the nilpotent homological operator for a 
nilpotent matrix iV is not as easy. There are two commonly used methods: 

• Given any scalar product (, ) on the space of polynomial vector fields, then an orthogonal 
complement to the range is coker Cn = kei^C^)*, the kernel of the adjoint. A nice 
choice of scalar product, generalizing the Frobenius inner product of matrices, leads to 
{C N f = C N * |Bel02llETB^87) . This is called inner-product style in [Mur03| . 

• From the representation theory of s((2) and the Jacobson-Morozov embedding theorem, 
for any nilpotent N there are matrices M and T that together give a representation of 
st(2). It then follows that ker Cm is a complement to the range of Cn [CS86llCS90b] . This 
is called st(2) style in |Mur03j . 

In general these approaches lead to different complements. Unfortunately, neither of these 
complements typically has the "simplest" form. For example, one could declare the form with 
the minimal number of nonlinear terms at each degree to be simplest. Murdock as discussed a 
separate procedure that can be appended to either style to simplify the normal form in this way 
|Mur03| . 

In this paper we study an incompressible vector field in three dimensions with a multiplicity- 
three eigenvalue A = so that J = + N. Generically A has geometric multiplicity one, so N 
has a single Jordan block. We will show that the complement to rng£ can be selected so that 
the normal form is particularly simple: 

Theorem 1. Consider a smooth vector field ([3]) on ]R 3 , £ = (x,y,z) T S M 3 , with vanishing 
divergence, (pQ) and linear part 



Dv(0) =JV= 1. (9) 




The vector field V can be transformed by a volume-preserving, near-identity transformation §5$) 
into the form 

N£ + e 3P {x,y), e 3 = (0,0,l) T . (10) 

up to terms of arbitrary degree. 

The simple complement proposed here extends Murdock's results to the case of divergence- 
free vector fields. For the particular N we are considering, the general result was already 
obtained in |ETB + 87j . It states that the nonlinear part of the transformed vector field has the 
form es(zipi+yip2+^P3) where tpi = (fi(x, y 2 — 2xz) are arbitrary polynomials. We will reconstruct 
this result in £j3j since we need an intermediate form from this construction to obtain Th. [TJ 
Indeed, as we explain in it is not possible to simply impose zero divergence on this simplified 
form. In other words, the operations of imposing zero divergence and computing the simplified 
normal form do not commute. 

A similar theorem holds for dimension two, i.e., the divergence- free case of the Takens- 
Bogdanov bifurcation. Generalization to four or more dimensions is not as straightforward, as 
we discuss in £J71 
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Normal forms for maps can be found by analogous means. Let / : M. n — > M n be a smooth 
diffeomorphism 

/(£) = J(£ + 6(£)) (11) 

where 6 = 0(2) represents the nonlinear terms. A formal conjugacy §5§ f = goip to & normal 
form g{rj) = J{r)+c{rj)) can be found if we can solve the homological equation Ljh(£) = c(£)— &(£) 
at each degree for h, where 

Lj = Adj-i -id (12) 

and Adjh(^) = J/i(J _1 £). The terms c in the normal form should again be selected to be in a 
complement to the range of the homological operator. Construction of this complement depends 
upon the properties of J. 

Here we consider the case of volume-preserving maps that have a fixed point with multipliers 
A = 1 of multiplicity three. The mapping analogue of Th. [1] is 



Theorem 2. Consider a smooth, volume-preserving diffeomorphism ()TTj) on E with the linear 
part 

fl 1 0\ 

Df(0) = J= 1 1 =I + iV. (13) 

\o o 1/ 

By a volume-preserving, near-identity transformation the map f can be put into the form 



g{Q = J{S + e 3 p(x,y)), e 3 = (0,0,l) T . (14) 

up to terms of arbitrary high degree. The truncation of the normal form at any degree is exactly 
volume preserving with a polynomial inverse of the same degree. 

The essential step for the transition from vector fields to maps is to use the exponential of 
the vector field and the exponential of ad instead of Ad. However, since the simplified normal 
form style is basis dependent, additional coordinate transformations are needed to transfer the 
result. We are motivated by the approach of [BC93j. 

These results can be extended to the (formal) unfolding of these bifurcations, as stated in 

Theorem 3. The unfoldings of the vector field of theorem 1 and the map of theorem 2 are 
obtained by replacing the polynomial p(x, y) in (jlOl) and (|14() by 

e + mx + fi 2 y + p(x, y; e, m, /i 2 ) , 

where the lowest order terms of p in x and y are quadratic. 



2 Polynomial Vector Fields 

Normal form theory uses formal power series of vector fields. In this section we introduce our 
notation for these vector fields and give some of the results about the polynomial subspaces and 
their bases that we will use. 

For each point £ € R n and each m S N n , let £ m = S,™ 1 !;™ 2 ■ ■ ■ £Jf n denote a scalar monomial 
with degree \m\ = Yli=l m i- Since the homological operators £ (|7|) and L (fT2j) preserve the 
degree, one can restrict to the finite-dimensional spaces of polynomial vector fields with fixed 
degree. 
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Let P n be the space of formal (without regard to convergence) polynomials in n-variables 
and C P n be those with homogeneous degree d, i.e., if p € then p(a£) = a d p(^) for 
all a G R. One basis for P£ is the set of all monomials {£ m : \m\ = d}. Consequently the 



dimension of PJ 1 is 



dimP d n 



n + d-l\ (n + d-l)\ 



d J " (n-l)!d! ' 

The vector space of polynomial vector fields in lR n is denoted F n and the subspace with entries 
from P'2 in each component is denoted 

F2 = {h(0:h i (OeP2,i = l,...,n}- (15) 
A basis for this space is given by all "vector monomials," 

K, m = re 4 = (0,...,0,r,0,...,0) T , \m\=d,l<i<n, (16) 

consequently, 

dimP d n = ndimP d ™. 

For example, dimPj = d + 1 and dimPj = 3(d + 2)(d + l)/2. It is usually obvious that the 
operators we consider that act on F n have F£ as invariant subspaces, i.e., they preserve the 
degree. 

The space of divergence free vector fields of degree-d is denoted by V£: 

V? = {v€F2 : V-v = 0}. (17) 

Since V- is a map from F£ to P&-\ the condition of vanishing divergence lowers the dimension 
by dimP^, and hence 

dimV^ = dimP£-dimP£_ 1 . (18) 
A basis for V£ can t> e constructed as follows: 

Lemma 4. A basis of V£ consists of n dim P^~ l basis vectors Pi^m where the " indicates that 
fhi = 0, \m\ = d, and (n — 1) dimP^_ 1 basis vectors of the form 

Vi, m = C 1 [(1 + m i+ i)iiei - (1 + mi)£ i+1 e i+1 ] , \m\ = d - 1 , (19) 

with 1 < i < n — 1 . 

Proof. First note that the dimensions of the two subsets add up: 

dim V£ = dim F r d l - dim P d n _ x = n dim P^ 1 + (n - 1) dim . 
The first set of basis vectors has vanishing divergence since 

vanishes when m% = 0. There are radimP^ -1 vector monomials of this type, since in the ith 
component of the vector we can put a polynomial depending on all n — 1 variables except 
and the dimension of this space is dimP™- 1 . The remaining basis elements are vector fields of 
the form A£ multiplied by £ m with \m\ = d — 1 and the diagonal matrix A = diag(ai, . . . , a n ). 
They are divergence free if 

n n n 

v • (cm) = E wcoiti) = E a ttM m + o = rE^K + 1 ) = °- 

i=l i=l i=l 
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For given m this equation has n — 1 independent solutions, those vectors that are orthogonal to 
the vector (1 + 77^)1=1,..,^. Hence there are (n — 1) dimP ( ^_ 1 such vector fields. The independent 
solutions can be chosen such that only two are nonzero, as in (|19|) . The A£ term gives 
monomials that depend on £j in component i, and are thus independent of the first group; they 
are also independent of each other because the monomials £ m are distinct. □ 

There is a natural inner product on the space of polynomial vector fields that generalizes 
the Frobenius inner product |Bel021lETB+87] . For each p,q€P2, define 

{p,q)=p(d^-q(Ok=o, (20) 

where "•" denotes the Euclidean scalar product and p(d^) is the polynomial p with each occur- 
rence of £j replaced by the derivative For example, for two vector monomials (116p . 

(Pi,mj Pj,rh) = 1Tl^Sm,rh^i,j j 

where ml = mi^! . . . m n l. It is easy to see that (p, q) = (q,p), and (p,p) > when p ^ 0. For 
linear vector fields, d = 1, this inner product reduces to 



(A^,B0=J2A ij B ij , 

the inner product that defines the Frobenius norm for matrices. 

In this norm, the two subsets of basis vectors for VJ 1 are orthogonal: 

Lemma 5. The basis vectors Pi,m, Vj,m °fVd are orthogonal for any i,j G [1, n] and m,fh £ N rf . 

Proof. This follows because pi^ is missing £j in the i th component, but Vj : m always has £j in 



component i if that component is non- vanishing. 
We define one additional subspace of F£: 

u n d = mm ■■ 0(0 6 pn ■ 

A basis for U2 is given by the set vector fields 



□ 



(21) 



re, 



i=l 



For example, 



771 



d-1. 





r 




(xy\ 


( xz 


C/f = span < 




xy , 




yz 






\xz ) 


\yz) 


U 2 



(22) 



Note that each basis vector in U2 has non-vanishing divergence V • u m = (n + d — l)£ m ; indeed, 
this space is the complement of the divergence free space. 

Lemma 6. For the inner product (j20|) . the spaces and V£ are orthogonal complements in 

rpn 

r d ■ 
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Proof. For any v E V£ and u = £ t/^ we have 

<e(Oe,«> = <fl,V-t;> = 0. 

Thus U2 and V^J 1 are orthogonal. That they are complementary spaces simply follows from the 
observation that their dimensions add to that of F£. By (|2ip . the dimension of U2 is that same 
as that of degree d — 1 polynomials, 

dimt/J = dimP^x . 

Therefore, (fT8|) gives dimF^ 1 = dim VJ 1 + dim[/J. □ 

Remark 1. For n = 3, the orthogonal decomposition = C/J seems reminiscent of the 
Helmholtz decomposition of vector fields, h = V x -0 + V^. However, it is different — in particular 
the curl of a vector 8£ in t/J does not vanish in general. The difference is that for the Helmholtz 
decomposition, the scalar product between vector fields is defined by the integral of the Euclidean 
scalar product (in M 3 ). Nevertheless, the dimension of the space of gradient vector fields that 
are not harmonic vector fields is dimPj +1 — (2(d + 1) + 1) = dimPj^, the same dimension as 

3 Normal form for vector fields 

In this section we will use the inner product style to compute the complement to rng£0 This 
will give the formal normal form for a vector field with a triple-zero eigenvalue. The results here 
reproduce those of [ETB + 87 . We will use an intermediate form of this result in the next section 



to project onto the divergence-free case and complete the proof of Th. [TJ 

The main advantage of the inner product (|20p is that it allows a simple construction of the 
adjoint of Cm = ad at on the space of polynomial vector fields |Bel02l lETB + 87 . Along the way, 



we prove the analogous result for the operator Adj that will be used for the map case. 
Lemma 7. Using the inner product (|20p . (Adj)* = Adj* and (ad^)* = ad^v*- 
Proof. From the definition (|20p 

(p,Adjq)=p(dt)-Jq(J- 1 Ok=o- 

Change coordinates to r\ = J -1 ^, noting that = J~ l *d v to obtain 

(p,Adjq) = (r P (J- u d v )) ■ q( V )\ v=0 = q(d v ) ■ (rp{J- l *n))\ v=Q = (g, Adj* p) . (23) 

The analogous statement for ad^ follows from Ad e tjv = e tadjv and (e tN )* = e tN * . Setting 
J = e tN in (|23|) and differentiating with respect to t at t = gives the result. □ 

This lemma implies that coker ad^v = ker ad at* is the orthogonal complement to the range of 
Cn- Using (jHJ), the cokernel of ad^v is therefore determined by the solutions of 

adAr. h = V N *h- N*h = , (24) 

where we introduce the linear operator T>n = N!^ ■ V. 



1 In our case this gives the same result as the s[(2) style 
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Here we solve (|24p for the n-dimensional generalization of Q 



N 



(25) 



that is, Ni j+i 



/o 1 ...\ 
1 ... 
1 

V : ■ o •■./ 

1, and Nij = otherwise. For the adjoint of N, the linear operator T> becomes 



V 



N* 



n-1 

i=i 
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i+1 



and 



is equivalent to the system of linear PDEs 



V N *hi = , V N *hj = Hj-i , j = 2, 3, . . . n . 



(26) 



This system is easily solved by the method of characteristics. The solution of the characteristic 
equations £ = N*£ and h = N*h is formally easy: 



h(t) 



e tN *h° 



(27) 



A solution of the PDE is obtained by inverting the equations for £(i) to solve for the invariants 



-tN* 



£. However, we must use one of the equations to eliminate time. The easiest way 



to do this is to use the equation for the second variable since it is linear in t. Denote the first 
three coordinates by (x, y, z) so that £ = (x, y, z, . . . , £ n ). Then ([2?]) gives y(t) = y° + tx°. When 
x(t) = x° ^ the characteristics all pass through the surface y° = 0. So up to this singularity 
we can choose this surface to define the initial conditions. To do this, set t = y/x, to obtain the 
invariants 



f(x,y,...U) 



(28) 



These invariants are rational functions with denominator x k ~ 1 . The first invariant is simply x 
itself, and the solution of the PDE will depend on arbitrary functions of the invariants. Thus we 
can clear the denominators in the remaining equations and define n — 1 polynomial invariants 



a 

P 
7 



2^3° = 

2,x 2 a 



2zx — y' L 
= 2>x 2 w — 



(29) 



3xyz + y , 

and so forth. Therefore, a formal solution to (126j) is 

h(0 = e y/* N *cj>(x,(3 n ,...), (30) 
of the n — 1 invariants. For example, for n = 2 the 



for an arbitrary vector valued function 
solution is 

h = 



<pi(x) 
f(/>i(x) +</> 2 (x) 



(31) 
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and when n = 3 we obtain 

^cP 1 ( x ,P) + ^ 2 (x,P) + ( f> 3 (x,P)J 



(32) 



The desired solution for h is in F?. It is, however, a nontrivial problem to obtain the most 
general polynomial solution. Even if we assume that the functions (pi are polynomials, the 
solution (|30p is generally rational because it contains powers of t = y/x. For n = 2 it is clear 
that 4>\ (x) must be x times a polynomial and the solution is simply 

h-( T ) (33) 

where cpi are polynomials that depend upon x only. 

For n = 3 one might a first think that (j>\ = x 2 <~p\ is required, since it appears with denom- 
inator x 2 in the third component of (|32|) . However, this is not the most general polynomial 
solution, since we can replace y 2 by 2zx — j3 to eliminate one power of x in the denominator and 
then remove the second by setting 0i = x<pi and <fo = <^y?i + </?3- This changes the third term 
to z(p% + ^1^2 + f3- Finally setting fa = x<f>2, we obtain the polynomial solution for n = 3 



h 



I Xtfil 

yipi + xip 2 ! • (34) 
\z(p! + y^ 2 + 9^3/ 



where (fi = (fi(x,(3) are polynomials of appropriate degree. This solution was shown to be the 
general polynomial solution in |ETB + 87] . 

With these forms it is easy to obtain the dimensions of the kernel of ad^* for two and three 
dimensions, giving the entries in the n = 2 and 3 rows of the leftmost pane of Table [TJ 

Lemma 8. The dimension o/kerad7v* in Pj is 2 and in P? is \3d/2 + 1]. 

Proof. This is a simple counting argument, based on the polynomial forms. For n = 2 at 
degvee-d, h contains two monomials, <~p\ = ax^ 1 and ip2 = bx d so the dimension is always two. 

For n = 3, the polynomials (fi in fl34j) whose arguments are of degree one and two, respec- 
tively, must be chosen appropriately. The even and odd cases can be treated separately: a 
degree-2m polynomial has the form 

m 

Y,a lX 2l (3 m - l ePL, (35) 

1=0 

and a degree 2m + 1 polynomial has the form 



m 

hx 2l+1 r~ l E P| m+1 • (36) 

1=0 



Each of these sums has m + 1 arbitrary coefficients, so the degree-ti case has [^rpl coefficients. 

To apply this to the form (f34"j) . note that when h is degree d, then 991 and c/?2 have degree 
d — 1. Thus each has [|] coefficients. The function (^3 is degree d, so it has [^tp] coefficients. 
Thus total number of coefficients is 
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Table 1: Dimensions of the various (co)kernels for the homological operator adjy with N in (|25p 
in dimension n at degree d. Th. Q]for n = 2, 3 implies that the first two rows of the middle and 
right tables coincide for any d. 

dim(keradjv) dim(ker adjv flV^ 1 ) dimP^ -1 



n\d 


2 


3 


4 


5 


6 


n\d 


2 


3 


4 


5 


6 


n\d 


2 


3 


4 


5 


6 


2 


2 


2 


2 


2 


2 


2 


1 


1 


1 


1 


1 


2 


1 


1 


1 


1 


1 


3 


4 


6 


7 


9 


10 


3 


3 


4 


5 


6 


7 


3 


3 


4 


5 


6 


7 


4 


7 


12 


17 


24 


31 


4 


6 


10 


14 


19 


25 


4 


6 


10 


15 


21 


28 


5 


11 


21 


36 






5 


10 


18 


31 






5 


10 


20 


35 






6 


15 










6 


14 










6 


15 











Since the each monomial in each function represents an independent vector in h, this is the same 
as the dimension. 

□ 

Remark 2. The construction of a general polynomial solution to (|26[) when n > 3 is complicated 
by the fact that there are polynomial combinations of the invariants that have a = x as a factor. 
For example the combination /? 3 + 7 2 = a 2 5 where 5 is a quartic polynomial in (x, y, z, w). Con- 
sequently, the power of x in the denominator of terms involving these combinations is different 
than its nominal value. It can be shown that is the only nontrivial relation in the case n = 4, 
and this leads to the polynomial form for this case [Mur03j . We are not aware that n = 5 has 
been explicitly worked out. Thus the dimension of these subspaces is harder to compute when 
n > 3 (see Table [TJ for n > 4). 

4 Divergence free vector fields 

In this section we complete the proof of Th. [TJ The remaining task is to compute a complement 
to the range of ad at for divergence-free vector fields. We start by showing that for any matrix 
N, the splitting F£ = V£ © 17% block diagonalizes ad^, that is 

ad N V2 C V? , and adjv 11% C 11% ■ (37) 

Moreover, since (adv)* = ad^v* this also applies to ad^y*. 

This block diagonalization allows us to compute cokeradAr for divergence- free vector fields 
by simply projecting the full coker onto V£ . This projection can be accomplished by solving 
an ODE for one of the components of h. However, the result is not the simplified normal form 
of Th. [TJ The last step is to show that the simplified form is also a complement to rngadAr by 
showing that its projection onto cokeradAr is of full rank. 

We begin by verifying (|37p in the next two lemmas. 

Lemma 9. The volume preserving subspace V£ is an invariant subspace of ad n for any matrix 
N; in other words, ad at V£ C V£. 
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Proof. Our goal is to show that if V-u = then V- (ad^r v) = 0. Writing the latter in components 
gives 

di{ad N v)i = di{N km £ m d k Vi - N ik v k ) = [Dv,N] u + (N£ • V)^ , 

where [ , ] is the matrix commutator and we use the summation convention. The first term is 
simply the trace of the commutator of Dv and N, but tr[A, B] = for any two matrices. The 
last term is (iV£ ■ V)(V • v), which vanishes since v € V£. □ 

Remark 3. Notice that it was not necessary to assume triV = in this Lemma because the 
trace of the commutator of any two matrices vanishes. For other Lie algebras, e.g. sp or so, it is 
necessary to assume that N is in the algebra to get a block diagonalization since otherwise the 
commutator [Dv, N] is not in the algebraH 

This lemma implies that the operator ad at acting on a vector in components v + u is block 
upper-triangular. To show that it is block diagonal, we must also show that ad at leaves the 
orthogonal complement {7? invariant. 

Lemma 10. The non-volume-preserving subspace U2 is an invariant subspace o/adjv for any 
matrix N. Thus ad N U2 C C/J. 

Proof. From (|2ip . the general element of U2 has the form where 9 £ ?2 is a scalar 

polynomial. Such a vector remains in U2 since 

ad^v 9£ = D(6£)N£ - N6£ = (V0 • N£)£ , 

is again of the form of scalar, degree-d polynomial times £. □ 

Thus coker ad^ can be split into two subspaces 

coker ad at = (coker ad^ (1U2) ® (coker ad at HVJ 1 ) , 

that are invariant under ad at. Thus, for a divergence- free vector field ([3]), the solution h to the 
homological equation ([6]) can be taken to be itself divergence free, and the normal form can be 
selected to be in coker ad^ nV d ™. 

We will only construct these subspaces for the two and three-dimensional cases. 

The two-dimensional case is simple. The general solution of the homological equation in this 
case was given in ([3T]) . Imposing V • h = implies that 4>i(x) = 0, so that h = e2<t>2(x)- This 
implies that dim coker ad at HV^ = 1, as shown in Table [TJ Note that the resulting vector field 
already has the "simplified" form of Th. [TJ In this case, we could also have imposed the volume- 
preserving condition on the polynomial solution f)33|) and obtained the same result. Neither of 
these two statements are true for the three-dimensional case 

Lemma 11. The intersection (coker ad^) D is the set of vector fields of the form 

h = xyipi + xij)2 j (38) 

\\y 2 ^i + y^2 + 1P3J 



2 Preserving symmetry in the normal form is different. There a matrix group G exists so that [G,N] = 0. 
This is replacing the condition on the differential of K and on N. When G is unitary this implies [G, N*] — 0. 
However, there is no associated invariant subspace, instead we find [Adc,adjv] = and [Adcadjv*] = 0, which 
means that the cokernel is invariant under the action of G. This implies that the normal form has the same 
symmetry group. However, this does not mean that there is a subspace invariant under adjv- 
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where ipi = ipi(x,[3) G Pj and ipz satisfies 



- 2dp^ 3 = xd x ^ + 3^i + Pdpih • (39) 

The dimension of the cokernel of adjv restricted to is d + 1 as in Table [7J 

Proof. An explicit computation of the divergence of (the not necessarily polynomial solution) 
([MP gives 

For any given <p\ and 02, this ODE for the dependence 03 on j3 can be integrated with respect 
to j3. For ^3 to be polynomial we require 0i = x 2 ip\. Similarly we need to set 02 = xip2 in order 
to make the 2 nd component of h polynomial. 

The polynomial ipi has degree d — 2, so by ([35]) and ([36]) it has pi^] coefficients. Similarly 
"02 is of degree d — 1 and has [~|] coefficients. For given and V'2 ) "03 is determined up to an 
arbitrary function of x, which introduces a single term cx d at fixed degree d. Thus the total 
number of coefficients is 

r d — 1, r cL 

r ^ _1 + r 2 1+1 = (i+1 - 

□ 

For completeness, we also can obtain a representation for the non-volume-preserving block. 

Lemma 12. The intersection (coker adjv)nt7j is the set of vector fields : 9(x,(3) £ -Pj_ 1 ,/3 = 
2zx — y 2 }. This space has dimension [|]. 

Proof. From Lem. [10] we see that we need to solve the equation 

ad^v* 9£ = (V0 ■ N*0 = V N «9 = ; 



which is the same as the first PDE in the system (I26p . The method of characteristics, as in §3] 
implies that the general solution is an arbitrary function of the invariants (|28|) as before. When 
n = 3 the invariants x and j3 have no combinations that have X clS Si factor, so the general 
polynomial solution is 

u = O(x,0)Z, 

for a polynomial function 6. Since 9 has degree d - 1, and (J3BJ, imply that it has [f] 
coefficients. □ 

Lemma [TT1 shows that that the dimension of ker ad^y HV^ is the same as that of the space of 
a single polynomial in two variables (given in the right pane of Tabled]). This makes it plausible 
that the divergence-free vector fields of the form h = (0, 0,p(x, y)) T form a possible complement 
of the range of ad^. This is proved in the next lemma. 

The simplified normal form was introduced by Murdock for the general case [Mur03j . He 
argued that a simplified complement to the range of adjy has the form of (|34p but with h\ = 
hi = 0. We show next that a similar projection can be done in the divergence-free case. 

Lemma 13. The set {e-sp(x, y) : p £ Pj} is a complement to rngad^ HVj. 



13 



Proof. We need to demonstrate that the orthogonal projection of the new set to the volume- 
preserving subspace coker ad at HV^ 1 is 1-to-l. In other words, the matrix of inner products of 
bases for the two spaces is nonsingular. 

Consider bases for the two subspaces. A basis for the simplified complement are the d + 1 
vector monomials {e 3 x k y d ~ k } where k = 0, 1, ... d. 

The basis for coker ad jv nV^ 3 given in (|38|) can be constructed by using monomials for the 
functions ip\ and rp2. For example when d = 2m is even, the monomial basis elements correspond 
to ipi = x 2 i~ 2 f3 m ~i and ip2 = x 2 - ?_1 /3 m_ - ? , for j = 1, 2, . . . m. The final basis element corresponds 

The d + 1 basis vectors are then 



to the free function obtained from integrating (139 j) . ifts 

( x 2j (3 m - j \ ( 



„2m 



x 



2j-l yf3 m-j 



V 2 



y 



m+j+l 
m—j+1 ' 



2j-2om-j 



h 



2,j 





x 2jprn-j 

\x 2 ^- 1 y(3 m ~^ 



The inner product of a general basis vector with h\ j gives 



I " ^k 2m— k t 

(e 3 x y ,h 



1,31 



i 

I k 2m— k i 

\x y ,^{y 



m + 1 



m 

k„,2m— 



m + j + 1 

JTT 



(5)x 



2j 







, h 3 = 


{"") 




\x 2m J 


- 2 (3 m ~i} 





m — j + 1 
2j — 2, that is when k 



k x 2j-2 



(-y 2 ) 



2\m-j+l\ 



This is clearly nonzero only for k = 2j — 2, that is when k = 0,2, 
functions ^2 j have inner products with the simplified basis given by 



, d — 2. The second set of 



(e 3 x k y 2m - k ,h 2 , j ) 



(x k y 2m 
(x k y 2m - k 



k x 2j-lfjm-j 



, X 



fin-]) 



which is nonzero only for k = 2j — l, or k = 1, 3, . . . d—1. The final inner product (e 3 x k y 2m ~ k , /13) 
is nonzero when k = 2m. Hence up to a permutation the matrix of scalar products between 
the basis vectors is diagonal and thus has full rank. Thus the orthogonal projection of the set 
e3p(x, y) to coker ad at is 1-to-l. A similar calculation pertains for the case that d = 2m + 1 
is odd; in this case the monomials in if>i cover the odd powers of x and the monomials in ip2 
cover the even powers. 

We have now proved Th. [TJ □ 

Remark 4. Note that while Murdock's simplified normal form (obtained by setting h\ = hi = 
in dMD) 

H X ,V, Z ) = h{z^pi{x,f3) +yip 2 {x,(3) +(p 3 (x,/3)) , (40) 

does form a complement to rngadjy for the general, non- volume preserving case, we cannot 
simply project the simplified normal form onto the volume-preserving subspace. Indeed, the 
only degree-d, volume-preserving elements for the simplified normal form are 



h(x,y,z) 



h(ay 2 x d - 1 + byx d - 1 + cx d ) 



giving a dimension of three, which is too small (except for d = 2). The point is that the 
transformation from inner product form to simplified normal form is equivalent to changing 
the scalar product. The new scalar product is defined by declaring the range of ad at and the 
proposed complement to be orthogonal. Assuming that the dimensions are right this is always 
possible. However, the new scalar product need not respect the divergence free conditions. More 
precisely, when the adjoint of ad at is defined with respect to this scalar product, then and 
U2 need not be invariant subspaces. 
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Remark 5. We have shown that ad at can be block-diagonalized into the divergence free subspace 
and its orthogonal complement. The basis given in £}2] achieves this block-diagonalization by 
Lem. EJ Thus it is interesting to compare the different representations of the complement of the 
kernel. We have just proved that another "simplified" complement is 

h = e 3 p(x,y) + 0(x,P)£ , 

which is different from (140p . Nevertheless, the projection argument works in this case as well. 

5 Volume preserving maps 

The normal form that we have just obtained for divergence- free vector fields can essentially be 
transferred to volume-preserving mappings by exponentiating ad at to obtain the homological 
operator for (|12[> . In this section we adapt the results of Bridges and Cushman [BC93] for 
Hamiltonian vector fields to divergence- free vector fields. The main difference is that in the 
volume-preserving case there is no scalar generating function. 

There are three basic steps to transferring Th. [T] to the mapping result Th. [2J First we note 
that the mapping homological operator L has the same adjoint properties as the flow operator 
C with respect to the inner product (f20j) . Then we use the fact that Ad is the exponential of 
ad to relate the fundamental spaces of these operators. Finally we show that it is possible to 
choose coordinates in which the linear form of the map is the exponential of the nilpotent block 
iV for the flow. 

As for the vector-field case, we will use the inner product style. Recall that Lem. [7] implies 
that (Adj)* = Adj* using the inner product (f20l) . Thus the mapping homological operator (fl~2l) 
satisfies 

(Lj)* = (Adj-! -I)* = Adj-i. -I = Lj* . 

Consequently, the construction of the normal form for the mapping (jlip reduces to finding 
a representation for kerLj*. Recall that Ad e jv = e adjv , thus there is a relation between the 
homological operators: 

Lj = Adj-i -I = exp(ad_ log j) - I , 

providing J has a (real) logarithm. This is true for (|14jl since none of the eigenvalues of J are 
negative. We will thus consider ad_i og j as the appropriate homological operator instead of the 
more difficult homological operator for maps. 

To use the relation between Ad and ad, we need to relate their fundamental spaces. The 
basic lemma relates the spaces for any nilpotent operator. This crucial step was inspired by 
[BC93] . 

Lemma 14. If N is a nilpotent matrix then M = exp N — I and N have the same kernel and 
range. 

Proof. Note that M = N(I + \N + . . . ^t^" 1-1 ) = NB = BN. B is a polynomial in N because 
N is nilpotent. We claim that B is nonsingular. Indeed when TV" is in Jordan form, then B is 
upper triangular and has a diagonal of all ones. Now, if w G rngM then there is a v such that 
w = Mv = NBv, so w G rng N. Moreover, ii w € rngiV then w = Nv = (NB)B~ 1 v = MB _1 v 
sow£ rngM. Thus rngM = rngiV. 

Similarly if z G kerM then = B~ x Mz = Nz and if z G keriV then = BNz = Mz. So 
kerM = keriV. □ 
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In order to apply the previous lemma to adjv we need the following well-known fact about 
adjoint representation of nilpotent matrices. 



Lemma 15. If N is nilpotent, then adjy is nilpotent. 

Proof. This follows from the fact that there exists a triad N, M, K that form a representation for 
sl(2), i.e., that [N,M] = K, [K, N] = 2N, and [K,M] = -2M and for any such representation, 
the matrices N and M are nilpotent and K is semisimple. Finally the operators ad/v, cloIm and 
adx also form such a representation. This proof is along the lines of Murdock [Mur03] Thm 
2.5.2, but it is a general theorem, see e.g. |Ser87| . □ 

The last step fixes the problem that Th.[T]is stated only for nilpotent Jordan blocks N, but 
log J = log (I + N) is not of this form. The remedy is to first do a linear change of coordinates 
that puts the volume-preserving map into a form where its linear part is exp N instead of 1+ N. 
Then Th. [2] connects to Th.[TJ The transformation must be chosen so that it does not destroy the 
simplified form of the map, in which all the nonlinearity is concentrated in the last component. 

Lemma 16. Suppose N is a nilpotent matrix, then there is an invertible matrix T such that 
T exp NT^ 1 = I + N. Moreover the conjugacy can be chosen such that T leaves e n invariant. 

Proof. While this lemma can be proved by induction for arbitrary dimensions, since we will 
apply it only for n = 2 and n = 3, we simply give the matrices for these cases. For n = 2 there 
is nothing to do since N% = 0, so that exp N2 = I + N%. For the 3x3 case we can see that 



exp N3 



I + N 3 + -Ni 



(I 


1 







1 


1 


Vo 








We have T 3 exp jV 3 = (I + A^Ts using 




□ 

Thus, using the new coordinate system, we can apply Th. Q] to the map case, giving the 
promised normal form (|14|) . 

To finish the proof of Th. [2] we only need to verify the claims about the properties of the 
normal form. Firstly every truncation of the map is exactly volume preserving. The determinant 
of the Jacobian of <?(£) = J(£ + e3p(x,y)) is the product of det J = 1 and det(I+ De3p(x,y)) = 1, 
since De~3p(x, y) is lower triangular. Another way of viewing this result is to say that time-one 
map of the vector field e%p{x,y) is the map £ + esp(x,y). Secondly the inverse of the map can 
be explicitly computed. Let £' = <?(£)• Then 



e3(z -y + z ,y - z) . 



This completes the proof of Th. [2l 
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6 Unfolding 



The unfolding of the triple eigenvalue collisions treated above is obtained by considering a family 
depending on e and additional unfolding parameters /Zj. The parameter e will be used to unfold 
the saddle-node bifurcation, while the parameters /ij will ensure that an arbitrary spectrum 
of the linearization can obtained near the bifurcation e = ^ = 0. The original map is now 
considered to depend upon e, and is expanded in a power series in (e,/z). The zero-order terms 
have been treated in the previous section. Normalization now proceeds order by order in (e, /i). 
The main observation (see, e.g. jETB + 87| ) is that at every order in the parameters the same 
homological operator is obtained. Thus all the results obtained above hold at each order, and 
hence instead of the polynomial p(x, y) with constant coefficients now each coefficient becomes 
a power series in (e, /x). 

The essential new feature is that now we also get a homological equation for terms of degree 
and 1 (previously these terms where absent by assumption). At degree the homological 
equation states that the constant vector must be in the kernel of the adjoint of N. The kernel of 
./V* is £3, which gives the constant term in the unfolding. At degree 1 the homological equation 
states that the linear terms must commute with the adjoint of N. The general matrix that 
commutes with ./V* is lower diagonal banded: 



Since we are dealing with divergence free vector fields we require /13 = 0. Now the second 
component of the vector field would have an entry Xfi?. To bring this into our simple form in 
which only the third component contains non-linear and unfolding terms we can simply add an 
element from the range of the homological operator (at degree 1). Thus we can modify this by 
adding any matrix (times £) that does not commute with N* . Hence we can simply remove 
e2Xfj>2, since it is in the range. In this way we find that the third component of the unfolded 
vector field becomes 

e + mx + fi 2 y + p{x, y; e, /ii, /i 2 ) . 

Under generic conditions on p one of the two parameters fj,i can be removed. Crossing over from 
divergence free vector fields to volume preserving maps works as before. This proves theorem 3. 

Interestingly the map thus obtained when p is restricted to be a quadratic polynomials is 
exactly the map that was studied earlier in [LM03J. In that paper this family of maps was 
studied because of its algebraic property that it is polynomial with polynomial inverse. We 
now see that this family of maps is also interesting because it is the unfolding of the volume 
preserving saddle-node bifurcation. 

7 Conclusion 

While we have obtained a normal form for three dimensions, our results do not apply to the 
higher dimensional case. Indeed, Table [H which compares the dimension of the kernel of the 
homological operator in the volume-preserving subspace and the dimension of the set spanned 
by £nP(£,i, ■ ■ ■ , Cn— l)j shows that these two dimensions are no longer the same when n > 3. It 
appears to be true that the dimension of the kernel is no larger and is strictly smaller for all 
degrees when n > 5. Thus we would conjecture that the span of e n p(£i, . . . , £ n -i) still contains 
a complement to the range so that that our normal form can still be used (though it would 





M2 
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not be the simplest normal form since it could contain more terms than strictly necessary). 
To prove this result the main problem is to obtain a description of the polynomial kernel of 
ad^-; this becomes very complicated when n > 3, see |Mur03j for the case n = 4. In higher 
dimensions most of the proof works, however, there is no simple general polynomial solution to 
the homological equation available as it was given in (|34[) . 

In a future paper we plan to study the dynamics of the unfolding of (|14p . 
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